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ABSTRACT 

We present the results of a program to acquire photometry for eighty-six late-M, L, and T dwarfs 
using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We examine the behavior of 
these cool dwarfs in various color-color and color-magnitude diagrams composed of near-IR and IRAC 
data. The T dwarfs exhibit the most distinctive positions in these diagrams. In M5.8 versus [5.8]- 
[8.0], the IRAC data for T dwarfs are not monotonic in either magnitude or color, giving the clearest 
indication yet that the T dwarfs are not a one parameter family in Toff . Because metallicity does not 
vary enough in the solar neighborhood to act as the second parameter, the most likely candidate then 
is gravity, which in turn translates to mass. Among objects with similar spectral type, the range of 
mass suggested by our sample is about a factor of five ('^70 MjupUer to '^15 M Jupiter), with the less 
massive objects making up the younger members of the sample. We also find the IRAC 4.5 /im fluxes 
to be lower than expected, from which we infer a stronger CO fundamental band at ~4.67 /xm . 
This suggests that equilibrium CH4/CO chemistry underestimates the abundance of CO in T dwarf 
atmospheres, confirming earlier results based on M-band observations from the ground. In combining 
IRAC photometry with near-IR JHK photometry and parallax data, we find the combination of Ks, 
IRAC 3.6 /im , and 4.5 /im bands to provide the best color-color discrimination for a wide range of 
M, L, and T dwarfs. Also noteworthy is the M^Ks versus ifs-[4.5] relation, which shows a smooth 
progression over spectral type and splits the M, L, and T types cleanly. 

Subject headings: infrared: stars - stars: fundamental properties - stars: late-type - stars: low-mass, 
brown dwarfs 



1. INTRODUCTION 

The study of low-mass dwarfs has progressed enor- 
mously in the last decade: from the strictly theoreti- 
cal beginnings for objects of sub-stellar mass (e.g. Ku- 
mar 1963; Grossman 1970; Nelson, Rappaport, & Joss 
1985) to the actual discovery in the late-1980's and mid- 
1990's of objects in two new spectral classes later than 
M (Beckhn & Zuckerman 1988; Nakajima et al. 1995). 
As a result of a large amount of observational effort over 
the past decade, the temperature sequence of dwarfs has 
been extended down to T ~ 700 K (Golimowski et al. 
2004a). Several hundred L dwarfs and several dozen T 
dwarfs have now been identified, leading to well-defined 
L and T spectral sequences (Kirkpatrick 2005), and to a 
large body of accurate photometry and quantitative spec- 
troscopy to characterize those spectral types (Leggett et 
al. 2001, 2002; Basri et al. 2000; Reid et al. 2002). Also 
of importance, other groups have been able to obtain 
accurate trigonometric parallaxes for a large number of 
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these very cool dwarfs (e.g., Dahn et al. 2002; Tinney 
et al. 2003; Vrba et al. 2004) and to obtain dynamical 
masses for a significant sample of stars (e.g., Henry & 
McCarthy 1993; Henry et al. 1999; Lane et al. 2001; Za- 
patero Osorio et al. 2004; Bony et al. 2004; Close et al. 

2005) . 

The observational progress was accompanied by an 
equally dramatic evolution in the sophistication of the 
theoretical models of cool dwarfs. The improved ob- 
servational data spurred theoretical model efforts such 
as Allard et al. (1994) and Brett (1995), who were 
able to incorporate much more extensive molecular line 
lists and thus improve the fit to observations and extend 
the model spectra further into the infrared. The model 
atmospheres have since rapidly evolved, incorporating 
ever-larger molecular line lists, dust formation and dust- 
settling, and improved treatment of pressure broadening, 
particularly for the alkali lines that dominate the optical 
and near-IR spectra of the L dwarfs. More recently, the 
models have begun to take into account non-equilibrium 
chemistry (e.g., Noll, Geballe, & Marley 1997; Saumon et 
al. 2000; Burrows, Marley, & Sharp 2000; Ackerman & 
Marley 2001; Allard et al. 2001; Lodders & Fegley 2002; 
Marley et al. 2002; Tsuji 2002; Saumon et al. 2003, 

2006) . 

While much of the work on sub-stellar mass objects 
to date has been carried out in the optical and near-IR, 
there is great potential to gain even more insight into 
brown dwarf atmospheres by extending observations into 
the mid-IR regime. The advantages include the peaking 
of the spectral energy distribution for L and T dwarfs 
between 1—4 /imand the presence of prominent molec- 
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ular features such as CH4, H2O, and NH3 (e.g. Bur- 
rows ct al. 1997; 2001). Ground-based infrared observa- 
tions are hampered longward of K band by telluric at- 
mospheric lines, such as OH and H2O, making the atmo- 
sphere increasingly opaque at these wavelengths. Ther- 
mal emission of the atmosphere and telescope structure 
also cause very high backgrounds for ground-based obser- 
vations longward of K band. Thus, space-based instru- 
ments with sensitivity in the mid-infrared are required 
for the study of sub-stellar mass objects in this wave- 
length regime. It was in this context that the Infrared Ar- 
ray Camera (IRAC) Guaranteed Time Observer (GTO) 
team decided that it would be important to obtain ac- 
curate mid-IR fluxes for a representative sample of M, 
L, and T dwarfs when the Spitzer Space Telescope was 
launched. 

We present in this paper the results of this survey. In 
section 2, we describe the sample selection of M, L, and T 
dwarfs and their basic properties. In section 3, we sum- 
marize the observing strategy used with Spitzer/IRAG 
and the subsequent data reduction and extraction of the 
IRAC photometry. In section 4, we present the basic 
color-color and color-magnitude diagrams for both IRAC 
and near-IR JHK photometry and discuss trending with 
color and spectral type. Finally, in section 5, we conclude 
with a few tentative interpretations of the observations. 

2. THE SAMPLE 

The targets in our program were selected from the lit- 
erature. Because the GTO program target selection was 
frozen a year prior to Spitzer launch on 25 August 2003, 
our selection was limited to those objects cataloged at 
that time. The primary consideration that governed our 
selection process was whether or not the target had a 
measured trigonometric parallax. The secondary con- 
straints on our selections were: (1) we wanted the sources 
to be relatively bright in order to ensure good signal-to- 
noise ratio on the photometry given our observing strat- 
egy (see section 3), since these objects would represent 
the fiducial sample for other science programs, (2) the 
targets needed to have well-determined spectral types be- 
cause of the desire to have representatives of all spectral 
sub- type bins from late-M through T types, (3) the tar- 
gets needed to be located in relatively uncrowded fields, 
and (4) should not be close binaries (separations of < 
6 arcseconds or about 5 IRAC pixels). While we tried 
to avoid including binaries in our sample, the informa- 
tion available in the literature was/is quite incomplete, 
so some now-known binaries were included in our pro- 
gram. Other objects may be binaries, but are currently 
not known as such. Given that dwarfs in the L7— T2 
spectral type range have been found to have a high bi- 
nary fraction (e.g. Burgasser et al. 2006), it seems in- 
evitable that more objects in our sample will turn out to 
be binaries than we have currently noted. 

The basic properties of the objects in our sample 
are summarized in Table 1. We provide here position, 
nomenclature, parallax (if available), spectral type, and 
near-IR photometry. Some target names arc abbrevi- 
ated from their official forms. Objects from the DE- 
NIS, 2MASS, and SDSS have had their survey acronyms 
shortened to three characters - DEN, 2MA, and SDS re- 
spectively - followed by the first four digits of right as- 
cension and then the first four digits of the declination. 



For the JHKs photometry, we use data from 2MASS 
primarily in order to maintain consistency. However in 
those cases where the 2MASS errors are large (> 0.20 
mag) or other flags in the database (e.g. blending, con- 
tamination, confusion) suggest questionable photometric 
reliability, we tabulate photometry from the literature 
(see references in Table 1). Because the near-IR pho- 
tometry for low-mass dwarfs is strongly system depen- 
dent (particularly for T dwarfs), we have chosen to rely 
on MKO-system photometry when reliable 2MASS pho- 
tometry is not available. We make this choice because 
there are well dociimcnted transformations between these 
two systems (Stephens & Leggctt 2004) as well as a 
wealth of MKO-system photometry for low-mass brown 
dwarfs available in the literature (see additional discus- 
sion in section 4.1.2). For two T dwarfs in our sample 
(SDS1624+0029 and 2MA1237+6526) we acquired new 
JHKs photometry with the PAIRITEL facility on Mt. 
Hopkins (telescope and instrument formerly used for the 
2MASS project, see Bloom et al. 2006 for additional de- 
tails). These photometry were calibrated to the 2MASS 
system using local comparison stars in each field-of-view. 
Because the T dwarfs arc intrinsically faint in the optical, 
it is difficult to classify these objects using the same crite- 
ria/features that work for the brighter L dwarfs. For the 
spectral types in Table 1, we have chosen to quote opti- 
cal types for the M and L dwarfs while, for the T dwarfs, 
infrared types are used. Near-IR classification systems 
take advantage of the peak of the spectral energy distri- 
bution being located in the infrared for T dwarfs and a 
wealth of molecular features in this wavelength regime. 
The unified near-IR classification system for T dwarfs 
(Burgasser et al. 2006) is compatible with optical clas- 
sification systems (e.g. Kirkpatrick et al. 1999, 2000) 
at type L8, creating a smooth continuum of types from 
M through T. It is worth noting, however, that for the 
L dwarfs, the optical and near-IR spectral types can be 
quite different, as they probe different optical depths in 
each wavelength regime (see discussion in Knapp et al. 
2004; Kirkpatrick 2005). 

In addition to our primary MLT sample, unsaturated 
photometry was secured for a number of late-M dwarfs 
from a related GTO program (PID 33 "A Search for 
Companions Around Stars Within Five Parsecs" ) . These 
include GJ 1002, LHS 288, GJ 412B, GJ 1111, LHS 292, 
SO0253-M652, and LP944-020. Archival data from an 
unpublished Spitzer Early Release Observation yielded 
photometry for the T dwarf binary e Indi BC, which are 
also present cxi in Tabk^ 1. Although GJ 229B was in- 
cluded in the nominal MLT GTO program, it was lost in 
the glare from GJ 229, so IRAC photometry could not 
be secured for this object. In total, supporting data for 
eighty-six objects for which we were able to secure IRAC 
photometry are presented in Table 1. 

2.1. Comments on specific sources 

As of the writing of this paper, the follow- 
ing objects have been determined to be binaries, 
unresolved or blended at the IRAC plate scale 
of ~1.2 arcseconds/pixel: GJ lOOlBC, SDS0423- 
0414, 2MA0746-F2000, GJ 337CD, SDS0926-F5847, 
SDS1021-0304, Kelu-1, 2MA1225-2739, 2MA1534-2952, 
2MA1553-F-1532, and 2MA1728+3948 (Reid et al. 2001; 
Burgasser et al. 2003a; Gizis et al. 2003; Bouy et al. 2004; 
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Golimowski et al. 2004b; Vrba et al. 2004; Burgasser, 
Kirkpatrick, & Lowrance 2005; Burgasser et al. 2005; 
Liu & Leggett 2005; Burgasser et al. 2006). The fol- 
lowing objects have been noted to have spectral pe- 
culiarities - unusually weak metal lines for their sub- 
type, emission lines in their spectra, or unusual col- 
ors - suggesting that some of these may be metal- 
poor subdwarfs, however, multiplicity may also play a 
role: 2MA0559-1404, 2MA0937+2931, 2MA1047-I-2124, 
SDSlllO+0116, 2MA1217-0311, 2MA1237+6526, and 
SDS1624+0029 (Burgasser et al. 2002; Burgasser et 
al. 2003b; Tinney et al. 2003; Gohmowski et al. 2004a). 

2MA0532+8246: This object, the first identified late 
type L subdwarf (Burgasser et al. 2003a), has been con- 
firmed as such in Reiners & Basri (2006). 

BRI0021-0214'- Photometric variability has been re- 
ported in /-band with a full amplitude about 0.05 mag 
and a possible periods of about 0.2 and 0.83 days (Martin 
et al. 2001). However, its large v sin i of ~34-42 km/s 
(Tinney & Reid 1998; Mohanty & Basri 2003), if true, 
would indicate a rotation period significantly less than 
0.8 days. 

DENIS0255-4700: Has a very high v sin i of 40 km/s 
(Mohanty & Basri 2003). A suggestion of photometric 
variability could be explained in part by rotational mod- 
ulation (Koen 2005). 

LP944-020: Listed as member of the Castor moving 
group by Ribas (2003), thus it's age is ~300-500 Myr. 
Tinney (1998) gives an age of 475-650 Myr for this star, 
based on the detection of lithium and the star's Ltoi- 

2MAO45 1-3402: Observed to have periodic photomet- 
ric variability with a period of '^6.9 hours (Koen 2004). 

2MA0746+2000AB: Reported to be photometrically 
variable with a period of 31 hrs according to Gelino 
(2002). Its measured v sin i (Bailer- Jones 2004) is in- 
consistent with that rotation period, however. 

2MA2224-0158: Noted as having an anomalously red 
spectrum, possibly indicative of an unusually thick con- 
densate cloud or of low surface gravity (Gushing et al. 
2005). 

3. OBSERVATIONS 

IRAG is a four-channel camera that obtains simul- 
taneous broadband images at 3.6, 4.5, 5.8, and 8.0 
/im (sometimes referred to as channels 1—4 respectively). 
Two nearly adjacent 5.2' x 5.2' fields of view in the fo- 
cal plane (for 256 x 256 pixel detector arrays) are viewed 
by the four channels in pairs (3.6 and 5.8 /im ; 4.5 and 8 
^m). Additional details on the design and performance 
of IRAG can be found in Fazio et al. (2004). 

Profiles of the IRAG filters are shown in Figure 1. For 
comparison, the profiles of the MKO L' and M' filters, 
as well as a model spectrum of a mid-T dwarf, are also 
shown to illustrate the relative band centers and band- 
widths of filters in recent use to characterize the infrared 
colors of M, L, and T dwarfs (e.g., Leggett et al. 2002; 
Golimowski et al. 2004a; Knapp et al. 2004). While the 
IRAC 3.6 /xm channel does share some similarities with 
the standard L filter, it is important to note that the 3.6 
and 4.5 //mlRAC channels are significantly bluer than 
the L', M, and M' filters and/or are much broader than 
all of these filters. The photometry presented in this pa- 
per have been calibrated using a zero magnitude defined 
by Vega in all IRAG channels (e.g. Cohen 2003; Cohen 



et al. 2003; Reach et al. 2005). As illustrated in Figure 
1, the spectra of low and sub-stellar mass objects have 
a great deal of structure within the filter bandpasses. 
Therefore, particularly for T dwarfs, IR fluxes measured 
using broadband fllters are critically dependent on both 
the central wavelength and the filter transmission profile. 
Significant color-terms can be present when transform- 
ing fluxes measured in apparently similar filters from two 
different systems (e.g., from L' to IRAC 3.6 /Ltm). 

3.1. Observing Strategy 

All observations in this GTO program share a com- 
mon Spitzer Astronomical Observation Request (AOR) 
design. Each target was observed using the same camera 
modes and frame times as well as dithering the targets 
over the same pixels (approximately) on the arrays. The 
primary motivation of this strategy was to minimize the 
introduction of additional uncertainties into the relative 
photometry for each object. Specifically, the AOR uses 
a 5-position Gaussian dither pattern, starting with the 
target near the center of the array. Dithering mitigates 
the effects of cosmic rays and bad/hot pixels, while the 
small scale factor option for the pattern keeps the tar- 
get near the center of the detectors in order to mini- 
mize spatially dependent uncertainties in the calibration 
of the instrument. The relative offsets for each dither po- 
sition in the pattern were identical for all of the objects 
in the program and were distributed within a radius of 38 
arcseconds of the center/initial position. High dynamic 
range mode was used because the program objects cover 
a wide range of brightness and the relative sensitivities of 
the four IRAG detectors preclude a single exposure time 
that will produce good signal-to-noise and unsaturated 
photometry in all four channels (see also discussion in 
section 3.3). The frame times of 30 and 2 seconds yield 
effective exposure times for each dither position of 26.8 
and 1.2 seconds respectively. 

The individual observations themselves were scheduled 
by the Spitzer Science Center (SSC) and were executed 
over the period 06 December 2003 to 01 November 2004. 
The data for each target, as received from the spacecraft, 
were processed by the IRAC pipeline software at the SSC. 
This pipeline removes the electronic bias, subtracts a 
dark sky image generated from observations of relatively 
empty sky near the ecliptic pole, flat-fields the data us- 
ing calibration observations of relatively blank fields near 
the ecliptic plane, and then linearizes the data using lab- 
oratory measurements of each pixel's response to a cali- 
bration lamp in frames of varying length^. The absolute 
calibration of IRAC is derived from aperture photome- 
try of standard stars. This calibration is applied to the 
data such that the final pipeline product for each frame 
is in units of surface flux per steradian (i.e.. Reach et al. 
2005). The pipeline processing produces calibrated data 
for each frame in the dither pattern for each IRAC chan- 
nel. This constitutes the Basic Calibrated Data (BCD) 
that were used for our analysis. 

3.2. Photometry 

Photometry was extracted for each source using the 
aperture photometry package in IRAF. To ensure the 

* Additional details can be found in the IRAC Data Handbook 
at http:/ /ssc. spitzer. caltech.edu/irac/data. html 
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highest signal-to-noisc and to minimize contamination 
by cosmic rays, the individual images in the dither pat- 
tern in each channel were combined using a local post- 
pipeline software suite developed by one of us (BMP). 
This software was used to co-register the frames in sky 
coordinates and then co-add the frames while rejecting 
temporally transient events (e.g., cosmic rays) as well as 
fixed-pattern noise and bad pixels. We note that caution 
must be used with any kind of filtering scheme to remove 
transients when co-adding or mosaicking IRAC data. Be- 
cause the point spread function (PSF) is under /critically- 
sampled by the IRAC arrays, this can lead to a pixel 
response function (PRF) that is sharply peaked if the 
PSF is centered on a pixel or to a PRF with a broader, 
lower-intensity peak if the PSF is centered at the corners 
of four adjacent pixels (the two extreme cases). Thus, 
the pixel phase (i.e., the position within a pixel) of the 
PSF centroid becomes important when using transient 
rejection based on, say, median-filtering and/or sigma 
clipping. The large changes in the peak of the PRF for 
a dithered source could lead to real flux being rejected 
and, thus, to an artificial reduction in the count rate of 
the source when combining the individual frames of the 
dither sequence (see additional discussion in Schuster, 
Marengo, & Patten 2006). 

To measure the photometry, we first removed the sur- 
face flux calibration from the images by dividing by the 
flux conversion factor (MJy/sr per DN/sec) found in the 
image headers of the BCD data and by multiplying by 
the efi'ective integration time - 26.8 or 1.2 seconds for 
the individual frames in the dither pattern and 134 or 6 
seconds for the co-added, cleaned data. The IRAC ab- 
solute calibration is based on observations of standard 
stars measured with aperture photometry using a source 
aperture with a radius of 10 native IRAC pixels in each 
channel. The background was estimated using an annu- 
lus centered on the source position with an inner radius 
10 pixels and width of 10 pixels (Reach et al. 2005). 
Because many of the targets in our program are located 
in semi-crowded fields, we chose to use a smaller source 
aperture with a radius of 4 native IRAC pixels in order 
to avoid contaminating flux from other, nearby sources. 
An additional benefit of using a smaller source aperture is 
an improved signal-to-noise ratio for many of the fainter 
sources. For the background estimation, we used the 
same annulus as that defined for the IRAC absolute cal- 
ibration. Photometry was extracted from both the co- 
added frames and the individual BCDs for each channel. 

The background-subtracted net source counts were 
then transformed into physical units by multiplying by 
the flux conversion factor and then by the solid angle 
for each pixel, to yield fluxes for each source. This flux 
density is then referenced against that of Vega in each 
IRAC channel to put our photometry on a Vega-relative 
system. Because we used a source aperture of 4 native 
IRAC pixels in radius, an aperture correction must be 
applied to the photometry. Because the correction for a 
4 pixel radius source aperture is not given in the IRAC 
Data Handbook, the aperture correction for the photom- 
etry was empirically determined by comparing our source 
aperture to the standard 10 pixel radius aperture using 
15 relatively bright, unsaturated targets from our pro- 
gram. Thus, the final calibration of our data used is that 
of Reach et al. (2005), combined with our own aperture 



correction. All of the relevant numbers used to calibrate 
our photometry are summarized in Table 2. Table 3 sum- 
marizes the IRAC photometry for all of our targets. 

3.3. Error Analysis 

The photometry for each of our targets were measured 
using both the co-registered, co-added, and cosmic ray 
cleaned mosaic images, and the individual BCDs in the 
dither pattern. We used the latter data to examine the 
repeatability of the photometry with a single AOR in 
order to gain some insight into errors associated with 
photometry of the co-added data. Figure 2 shows the 
standard deviation of the photometry from the individual 
BCDs for each target for each IRAC channel. All of our 
targets were observed using HDR mode, meaning both 
a short and long frametime exposure was taken at each 
position in the dither pattern. In general, we favored the 
use of the long exposure (30-second frametime) BCDs 
over the short exposure (2-second frametime) BCDs in 
our HDR AORs, unless the target object was saturated 
in the long exposure frames. In a few cases where there 
was good signal-to-noise in both the long and short ex- 
posure frames, we found good agreement (within a few 
percent) in the calibrated photometry from both frame- 
times. We find that the mean fluxes from the individual 
measures compare very well to those measured in the mo- 
saic images. In fact, we used this comparison as one of 
our checks of photometric consistency in detecting out- 
liers in photometry of our target in the individual frames, 
induced by bad pixels or cosmic ray hits. 

Overall, we find that the dispersions of the individ- 
ual measures used to create the mosaics are larger than 
those estimated using photon statistics and basic detec- 
tor characteristics (e.g., read noise, gain, etc.) alone. For 
IRAC channels 1 and 2 these are 5 and 3 times larger 
respectively, while for channels 3 and 4 they are essen- 
tially comparable, with the individual measure disper- 
sions being only ~20% larger on average in each channel 
than the basic photon statistics. This is probably due in 
part to some camera properties that have not been fully 
characterized in the current calibration. These include 
intra-pixel phase sensitivity variations, variations in the 
pixel solid angle due to geometric distortions introduced 
by the telescope and camera optics, and variations of 
the spectral response over the arrays due to the tilt of 
the filters with respect to the optical path and the rel- 
atively wide field of view of the camera. For the latter 
two issues, the corrections are defined to be unity at the 
center of the arrays (Reach et al. 2005). Thus, while all 
of these effects will act to increase the dispersion of the 
individual measurements in our dithered data, because 
our observing strategy for this program keeps the target 
object within ~38 arcseconds of the center of the arrays, 
these additional uncertainties should be minimized. 

We have chosen to use the standard deviation of the 
individual BCD measures for each source in each channel 
as the uncertainty of the photometry. These are the un- 
certainties quoted in Table 3. The number of BCDs used 
to construct the co-added frame for each channel and to 
calculate the standard deviation are also listed for each 
source and channel in Table 3. In terms of absolute cal- 
ibration to the Vega-relative magnitude system, the un- 
certainty of the IRAC photometry, convolved with the 
additional uncertainty introduced by our aperture cor- 
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rection, is about 3% in IRAC channels 1—3 and 5% in 
channel 4 (Reach et al. 2005). 

4. IRAC MAGNITUDES AND COLORS FOR LATE-M, L, 
AND T DWARFS 

Molecular features found in the infrared for low mass 
stars and sub-stellar mass objects have strengths that are 
a strong function of temperature and pressure. With the 
appropriate choice of photometric bandpasscs, one can 
study trending as a function of color, magnitude, and 
spectral type as a part of characterizing the atmospheres 
of these objects. For example, in the near-IR, the spec- 
tra of L dwarfs are characterized by absorption from CO 
and H2O, while the T dwarfs arc dominated by broad ab- 
sorption bands of CH4 and H2O, as well as collisionally 
induced absorption (CIA) by H2. Thus in near-IR pho- 
tometry, M and L dwarfs become redder with decreasing 
Teff in J — H and H — K. The L to T dwarf transition oc- 
curs as the silicate and iron condensates (clouds) become 
buried at increasing depth in the latc-L dwarfs. II2O ab- 
sorption begins to dominate the near-IR spectrum, lead- 
ing to a bluing of the near-IR colors through the early-T 
types. The colors then become even bluer from early-T 
to late-T with the onset and growth of CII4 absorption 
and CIA II2 ui K . The overall result is that the J — H 
and H — K colors for T dwarfs become bluer with in- 
creasing spectral subtype, becoming degenerate with the 
colors of higher mass K and M dwarfs. 

The IRAC filters were selected primarily to provide 
contiguous bandpasses from ~3 to 10 /xmfor the deter- 
mination of photometric red-shifts for extragalactic ob- 
jects as part of the primary Spitzer mission objectives 
(Fazio et al. 2004). However, these bandpasses were also 
defined, in part, to provide diagnostics for the study of 
sub-stellar mass objects. The bandpass for IRAC channel 
1 includes much of the CII4 fundamental absorption band 
(~3.3 /Lim). Channel 2 includes the continuum peak that 
is present for all stars cooler than 3000 K, making this 
the most sensitive IRAC channel for the study of sub- 
stellar objects. Channel 2 also contains the broad but 
shallow CO fundamental absorption band (^4.7 /im), 
whose presence in the T dwarfs provides evidence for 
non-equilibrium chemistry models (see discussion in sec- 
tion 5). Channel 3 includes II2O absorption. Finally, for 
channel 4 the most important molecular absorption in 
this bandpass is due to CII4. 

IRAC provides very precise fluxes in a wavelength 
regime that is poorly studied to date. For this reason, we 
believe it is useful to examine the IRAC data from several 
perspectives. First, we take a purely empirical viewpoint 
and limit ourselves to only the IRAC photometry to ex- 
plore the IRAC colors of M, L, and T dwarfs. Next, 
we combine IRAC and near-IR photometry to examine 
the location of M, L, and T dwarfs in both color-color 
and color-magnitude diagrams. In section 4.3 we show 
the variations of near- and mid-IR colors with spectral 
type. Finally, in section 5, we compare color-magnitude 
diagrams for M, L, and T dwarfs with theoretical models. 

4.1. Color-Magnitude and Color- Color Diagrams 

4.1.1. IRAC Photometry 

Star-formation studies often use plots of [3.6]-[4.5] ver- 
sus [4.5]-[5.8] or [3.6]-[4.5] versus [5.8]-[8.0] in order to 
identify young stars with warm circumstellar dust disks 



and envelopes. This technique is particularly effective be- 
cause normal stars have essentially Rayleigh- Jeans spec- 
tra at IRAC wavelengths, and hence have colors very 
near zero. Stars with warm dust disks have significantly 
red colors in IRAC wavelengths, and thus separate from 
older, dust-free stars. Late-M, L, and T dwarfs, how- 
ever, depart from this convenient scenario because they 
are cool enough to have photospheres sufficiently pol- 
luted by molecules that their spectra depart greatly from 
black-bodies. 

As shown in Figure 3, L dwarfs have [3.6]- [4.5] colors 
that remain close to zero, but have [4.5]- [5.8] colors that 
become significantly bluer. For the T dwarfs, the strong 
redward trend in [3.6]-[4.5] is the result of CH4 absorp- 
tion removing increasing amounts of flux from the 3.6 
/xm bandpass as Toff decreases. On the other hand, [4.5]- 
[5.8] trends strongly blueward for the T dwarfs, reflect- 
ing the strengthening of H2O absorption with decreasing 
Tcu in the 5.8 ^m bandpass, relative to the CO absorp- 
tion in the 4.5 /xm bandpass. 

For [3.6]-[4.5] versus [5.8]-[8.0] (Figure 4), the M and L 
dwarfs show about the same range of color in [5.8]- [8.0] 
as in [4.5]-[5.8], but in this case there appear to be more 
overlap in the dispersion between these two types. Again 
though, it is the T dwarfs that stand out in this plot, 
with both the [3.6]- [4.5] and [5.8]- [8.0] color indices trend- 
ing strongly redward for decreasing Teff . For the latter 
color index, the H2O feature in the 5.8 /xm bandpass ap- 
parently decreases the flux in this filter faster than the 
weaker CH4 absorption band in the broader IRAC 8.0 
^m bandpass. 

To further explore the trending in the IRAC colors, we 
provide several variants on one specific IRAC color-color 
diagram - a plot of [4.5]-[5.8] versus [3.6]-[8.0]. In the 
first of these plots. Figure 5a shows our sample (exclud- 
ing binaries and spectrally peculiar dwarfs) with sym- 
bols corresponding to the object's spectral type. Fig- 
ure 5a illustrates that the three spectral types inhabit 
relatively distinct parts of color-color space, with the 
[3.6]-[8.0] color becoming progressively redder for later 
spectral types (i.e., that it correlates with effective tem- 
perature), whereas the [4.5]-[5.8] color is approximately 
constant for M and L dwarfs, but turns sharply blueward 
with later T types. There is not a one-to-one correspon- 
dence between spectral type and IRAC color however - 
the M and L types overlap in their IRAC colors, whereas 
the T dwarfs appear to have a very large scatter in their 
colors. Are these effects simply "noise" , or are they in- 
dicative of complexity in the spectra of these stars be- 
yond that attributable to just spectral type? 

To address these questions, we construct several other 
versions of Figure 5a but now isolating specific spectral 
subtype ranges. Figure 5b shows the location of "early" 
and late M dwarfs in this color-color plane, and illus- 
trates that at least for our sample of stars, field M dwarfs 
with spectral type MS and earlier have IRAC colors es- 
sentially consistent with zero, with only the M8.5 and 
later M dwarfs having significantly non-zero IRAC col- 
ors. It is only the stars in the latter sub-type range that 
overlap in color with the L dwarfs. Figure 5c similarly 
isolates several L dwarf subtypes. We have combined the 
LO through L3.5 dwarfs into one group because we see 
no significant variation in IRAC color within this subtype 
range; however, the IRAC colors of these early L dwarfs 
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are significantly different from the colors of M dwarfs 
with type earlier than M8. The L4 and L4.5 dwarfs ap- 
pear to be a transition class, with two dwarfs of this type 
having IRAC colors like those or earlier type L dwarfs, 
while the other two are significantly redder, particularly 
in their [4.5]-[5.8] color. The L7.5 to L8 dwarfs have 
redder colors than earlier L dwarfs, and have [4. 5] -[5. 8] 
colors that are redder, on average, than even the early T 
dwarfs. 

Figure 5d highlights the IRAC colors of T dwarfs. The 
early T dwarfs have [3. 6] -[8.0] colors that overlap with the 
latest L dwarfs, but they are displaced blueward in [4.5]- 
[5.8] color from those dwarfs. In general, the T dwarfs 
show a progression, with their [3.6]-[8.0] color becom- 
ing redder and their [4.5]-[5.8] color becoming bluer with 
later T subtype. The T5 and T6 subtypes appear to be 
transitional, with a relatively large spread in [3.6]- [8.0] 
color for the T5 dwarfs and, in comparison, colors for 
the two T6 dwarfs in our sample that are closer to those 
of the T4 dwarfs than the T7-T8 dwarfs. 

Finally, in Figure 5e, we show the relative colors of 
those objects in our sample that are considered spec- 
trally peculiar or unusual. These include the lone sdL 
in our sample (2MA0532-h8246) and the lone T dwarf 
with a peculiar designation in its spectral type (T6p, 
2MA0937-I-2931). The other T dwarfs in the figure have 
been noted in the literature as being spectrally unusual. 
For the latter, how much of their peculiarity is due to 
atmospheric anomalies (i.e., metallicity) or multiplicity 
remains to be seen. For example, only recently has 
the unusual spectrum of 2MA0423-0414 been revealed 
to be a composite spectrum of an plus '^T2 close 
binary system (Burgasser et al. 2005). Nevertheless, 
these objects as a whole; seem to have discrepant colors 
when compared to other L and T dwarfs, particularly 
for 2MA0532-I-8246 (sdL), 2MA0937-I-2931 (T6p), and 
2MA1217-0311 (T7.5). 

4.1.2. IRAC Photometry Combined with JHK 
Photometry 

The combination of IRAC broadband photometry with 
near-IR photometry^ breaks the degeneracy of the M and 
L dwarfs colors. In a plot of iir5-[3.6] versus [3.6]- [4.5] 
(Figure 6), the iiTs- [3.6] color index shows a clear trend 
redward for the M and L dwarfs, well correlated with 
spectral type. For the T dwarfs, this same index trends 
over a more narrow range of color for increasing subtype. 

^ For this comparison to work, it is necessary for all of the near- 
IR JHK photometry to be on the same photometric system. While 
a number of JHK systems arc reasonably similar, with only a few 
percent difference in the fluxes for observations of normal stars 
(e.g., Leggett et al. 2002, Bessell & Brett 1988), these diflterences 
are magnified for low mass and sub-stellar mass objects, where 
strong molecular absorption bands in the infrared lead to JHK 
magnitudes that are highly dependent on the bandpasses of the fil- 
ters used. The recent introduction of the MKO near-IR photomet- 
ric system (Tokunaga et al. 2002) has made this problem especially 
acute for the L and T dwarfs, where differences of --^20% or more 
are observed between MKO photometry and that of other JHK 
systems for the same object. For all figures in this jjaper that use 
JHK photometry, we have transformed these data to the 2MASS 
system using relations from Stephens & Leggett (2004). While wc 
recognize that the MKO system has been endorsed by the lAU as 
the preferred photometric system for ground-based ncar-IR obser- 
vations, the vast majority of the JHK data for the objects in our 
sample are already in the 2MASS system. 



Although the i^5-[3.6] color index is almost degenerate, 
the [3.6]- [4.5] color index still spreads them nicely in Fig- 
ure 6. Because pressure induced H2 absorption affects 
the Ks band in particular and is stronger if gravity is 
high or metallicity is low, this can introduce scatter in 
Ks- This is likely the major contributor to the large 
scatter in color seen with the same spectral subclass in 
the mid- and late-T dwarfs (see also section 4.3 for the 
trending of color versus spectral type for these two color 
indices). For Ks-[4:.5] versus [3.6]-[4.5] (Figure 7), the 
same breaking of degeneracy and range of color is seen for 
the M and L dwarfs in if5'-[4.5] as for _fi'5-[3.6]; however, 
in this case the T dwarfs show a strong redward trend, 
with a spread of '~2.5 magnitudes in Ks-['i.5]. The rela- 
tive amount of scatter in the color index for the mid- and 
late-T dwarfs remains the same as that seen for 7^5- [3.6] 
in Figure 6. 

4.2. Absolute Magnitude versus Color 

For those objects in our program with trigonometric 
parallax measurements, we present two color-magnitude 
plots that could be potentially useful for identifying low- 
mass dwarfs in the field or as companions to other stars, 
using a combination of IRAC and near-IR photometry 
(some additional IRAC color-magnitude diagrams focus- 
ing on the T dwarfs are also presented in section 5). 

In Figure 8, versus i4'5-[4.5] shows that the colors 
of M, L, and T dwarfs follow a relatively smooth progres- 
sion with increasing type. This is similar to the trending 
seen for M^s versus Ks-M' as reported by Golimowski 
et al. 2004a. On the other hand, in Figure 9, Mks 
versus J- [4.5] shows almost 2.5 magnitudes of range for 
the M and L dwarfs, while for the T dwarfs the J-[4.5] 
color is essentially degenerate with the L dwarfs (and 
with themselves, the early-T dwarfs falling in the tran- 
sition area between the L dwarf and T dwarf plateaus, 
being degenerate not only with the mid-L dwarfs but also 
with mid-T types). It is only the intrinsic difference in 
luminosity that lifts the degeneracy for the J- [4.5] color 
index. 

4.3. Color versus Spectral Type 

In the same spirit as Leggett et al. (2002), Knapp et 
al. (2004), and Golimowski et al. (2004a), we summarize 
the trending of various color indices against spectral type 
(Figure 10) as a series of "postage stamp" plots. While 
most of this trending has already been detailed earlier in 
this section in the color-color and color-magnitude dia- 
grams, plotting the indices against spectral type reveals 
some common themes. In particular, most of the plots 
in Figure 10 show breaks in the trending near the L/T 
boundary. The most dramatic example is in [3.6]- [4.5] 
where the slow blueward trend through the M and early- 
L types turns redward at mid-L (given the appearance 
of CH4 absorption in the IRAC 3.6 /xm bandpass) and 
heads strongly redward through the T types. On the 
other hand, the _ftrs-[4.5] and [4.5]- [5.8] plots show breaks 
at the mid-T types. For [4.5]- [5.8], the blueward tmn at 
mid-T presumably happens due to the strengthening of 
H2O absorption with later spectral type in the IRAC 5.8 
^m filter, while for Ks-[4:.5], the onset of both CII4 ab- 
sorption (2.2 /xm overtone feature) and CIA H2 in the 
Ks band drives the color index redward for later types. 



IRAC Photometry of M, L, and T Dwarfs 



7 



The red ward trend of the _ftr5-[3.6] color index turns blue- 
ward as the fundamental CH4 absorption feature at 3.2 
/xm switches on at mid-L, and then turns redward again 
at mid-T when this feature saturates and 2.2 /xm CH4 ab- 
sorption and CIA H2 in the Ks band assert themselves 
at mid-T. 

These color versus spectral type plots serve to show 
that in the previously presented color-color diagrams 
based on IRAC photometry alone, much of the scatter 
in the T dwarf colors come from the [4.5]- [5.8] and [5.8]- 
[8.0] color indices. The trending of color versus spectral 
type for the [3.6]-[4.5] color index is very smooth. While 
some of the scatter in the [5.8]- [8.0] color index could be 
due to the interplay of H2O absorption (starting in the 
^late-M types) versus the 7.7 /im overtone band of CH4 
in the IRAC 8.0 /xm filter (beginning at ~mid-L types), 
it is also true that our fixed integration time AORs yield 
the highest quality data for IRAC channels 1 and 2, the 
most sensitive of the four IRAC channels, and thus sys- 
tematically lower signal-to-noise data for our sample in 
IRAC channels 3 and 4 (i.e., some of the scatter seen in 
[5.8]- [8.0] might be due to the larger errors on this color 
index compared to the others presented in Figure 10). 

5. DISCUSSION 

As shown in the previous sections, the infrared colors 
of low mass stars and brown dwarfs are generally well- 
correlated with their spectral types. In turn, the spec- 
tral types of M, L, and T dwarfs are reasonably well- 
correlated with their effective temperatures, though to a 
lesser degree than for higher mass objects because molec- 
ular chemistry and particulate clouds compete aggres- 
sively in determining the detailed spectral energy distri- 
butions of cool dwarfs. Theoretical models are beginning 
to be able to match the observed spectral properties of 
very cool dwarfs, and we have used those models to help 
interpret the IRAC photometry in Section 4. 

While a detailed theoretical analysis of these new 
Spitzer/'IKA.C data is beyond the scope of this paper, 
there are a few conclusions of a general nature that can 
be extracted at this preliminary stage. To do so, we 
have generated three color-magnitude plots (Figures 11- 
13) for most of the T dwarfs in our sample with par- 
allaxes and superposed synthetic magnitudes and colors 
derived from theoretical spectra for an effective temper- 
ature (Tefl ) range from 700 K to 1300 K in steps of 100 
K at three gravities. These gravities are 10^ ''' cm s~^ 
(green circles), 10^ cm s~^ (red triangles), and 10^'^ cm 
s~^ (blue squares). Figure 11 depicts the absolute mag- 
nitude in the IRAC 3.6 /um bandpass versus the [3.6]-[4.5] 
color. The legend on this figure indicates the T dwarfs 
included on this figure and they are numbered in order 
of spectral sub- type. 

As Figure 11 demonstrates, despite the fact that the 
numbering is in order of increasing spectral sub-type, the 
3.6 /xmlRAC data are not monotonic in either magnitude 
or color. That they are not in order in magnitude may 
in part be ascribed to multiplicity. However, this cannot 
be the; (explanation for the majority of these T dwarfs. 
Furthermore, the [3.6]-[4.5] colors too are non-monotonic 
in spectral type. This is even more obvious in Figure 
12. which portrays M5 g versus [5.8]-[8.0]. For example, 
SDS1346-0031 and 2MA0727-M710 only differ by half a 
spectral sub-type, but they have very different colors. 



The non-monotonicity of these Spitzer/l^AC data is 
the clearest indication yet that the T dwarfs are not a 
one-parameter family in Tgff , but that more than one 
parameter is influential in determining the spectroscopic 
type. Because the metallicity can not generically vary 
enough in the solar neighborhood to explain this, the ex- 
tra parameter may be gravity. This was the conclusion 
of Burrows et al. (2002) and Knapp et al. (2004) us- 
ing different datasets at shorter wavelengths, and is con- 
firmed here. For the T dwarfs, gravity will translate into 
mass (Burrows et al. 1997). From the comparison of the 
spread in the data in Figures 11 and 12 with the spread 
with gravity of the theoretical models, we conclude that a 
range of gravities of approximately a factor of five is rep- 
resented in the extant T dwarf sample. This can be con- 
verted into a range of masses from ~70 yijupiter to ~15 
^^.jupiter- The less massive objects would also be younger 
and a range of ages of about a factor of ten (0.2—0.3 to 
10 Gyr) is indicated (Burrows et al. 1997; Burrows et al 
2001). The latter conclusion might be problematic, but 
such is indicated by our preliminary analysis. 

Figure 13 is an HR diagram of IRAC M4.5 versus the 
IRAC [4.5]-[5.8] color. The non-monotonicity of the data 
seen in previous plots color survives in this plot as well. 
However, as is abundantly clear in the figure, the theo- 
retical models do not fit these data. Because CO has a 
strong spectral feature at ^4.67 /im , one can interpret 
the discrepancy between theory and the IRAC data as an 
indication that equilibrium CH4/CO chemistry underes- 
timates the abundance of CO in T dwarf atmospheres. 
This conclusion was already reached by Golimowski et al. 
(2004a) using M-band measurements from the ground. 
The magnitude of the discrepancy translates into an over- 
estimate in the ~4-5 /imflux by factors of 1.5 to 3.0. 

It should be noted that the opacities of water in 
the mid-infrared are still being studied and refined 
(Schwenke 2002). Continuing ambiguity in the H2O 
opacities are a source of systematic theoretical error, 
with the product of the O abundance and the H2O opac- 
ity directly connected to the goodness of fit at 5.8 and 
8.0 /zmfor the later T dwarfs. Moreover, at the low- 
est Teff represented here, the condensation of K into KCl 
and Na into Na2S could introduce hazes with interesting 
optical depths. Finally, the early T dwarf spectra show 
signs of silicate clouds (Stephens, Mar ley, & Noll 2001; 
Marley et al. 2002; Burrows et al. 2002) not included in 
the models we have presented. Hence, the study of sub- 
stellar mass object atmospheres is still in its formative 
years. 

Nevertheless, as Figure 14 implies, one can be encour- 
aged that the basics are falling into place. Figure 14 com- 
pares a theoretical T dwarf spectral model with T^g /g 
= 750 K/IO'"^ cm s^^ (Burrows et al. 2002; Geballe et 
al. 2001) with the four IRAC fluxes measured for the 
T7.5 dwarf GJ 570D. This model was generated in 2002 
for GJ 570D to fit its optical spectrum shortward of 1.0 
/xm (Burrows et al. 2002). A slight (~40%) discrepancy 
in the 4. 5-/^m bandpass, attributable to a CO abundance 
excess in its atmosphere, is visible. Despite this, Fig- 
ure 14, and Figure 7 in the related paper by Burrows et 
al. (2002), together represent an acceptable fit from 0.6 
jumto ~8.0 /umand indicate that the IRAC data were 
successfully anticipated. 
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Fig. 1. — A comparison of the IRAC bandpasses with the MKO L' and M' filters in recent use for ground-based studies of M, L, and T 
dwarfs. Also shown for comparison is a model spectrum for a Teg = 950K, g = 10^ cm s~-^ ~mid-T dwarf (Burrows et al. 2002) and that 
for an AO V star. The IRAC bandpasses have been scaled to show their relative sensitivity. 
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Fig. 2. — Repeatability of photometry for targets versus their cahbrated magnitudes for the four IRAC channels. Shown here are 
the standard deviations of the individual observations in the AOR dither pattern (five for most objects) versus the calibrated magnitude 
derived from the co-addition of those individual observations. For the 5.8 and 8.0 fim channels, the 3cr limiting magnitude for a high 
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Fig. 3. — [3.6]-[4.5] versus [4.5]-[5.8] color-color diagram for all of the sources in our M, L, and T dwarf sample, excluding known binaries 
and spectrally peculiar objects. As indicated in the legend, different plot symbols and colors are used to represent the M dwarfs {yellow 
diamonds), L dwarfs {green squares), and T dwarfs {red triangles) . The error bars to the right of each plot symbol in the legend represent 
the median values for each spectral type for the color indices used in this figure. 
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Fig. 5a. — [3.6]-[8.0] versus [4.5]-[5.8] color-color diagram, (a) Same plot symbols and objects as plotted in Figure 3. (b) Same data 
as in (a) except colored symbols show region of color-color space occupied by spectral types <M8 {yellow squares) and MS— M9 {green 
circles), (c) In this case the colored symbols show a clear separation in color-color space for early-L types (LO— L3) {blue squares) and 
late-L (L5— L8) {red triangles). The L4 dwarfs {cyan circles) appear to form a transition class between the early and late types, (d) The 
break out of the T dwarfs shows that for the early types TO— Tl {blue squares), T2— T3 {cyan triangles), and T4 {crosses), there is a trend 
towards redder [3.6]-[8.0] color with increasing spectral type, but little discernible change in the [4.5]-[5.8] color over these same types. The 
late-T types T7— T8 {red asterisks) are much redder in [3.6]-[8.0] and bluer in [4.5] — [5.8] than the early types. The T5 {green circles) and 
T6 {yellow diamonds) subtypes appear to be transitional, with a relatively large spread in [3.6]-[8.0] color for the T5 dwarfs and redder 
[4.5]-[5.8] colors for the T6 dwarfs, (e) Spectrally peculiar objects in our sample. Each object is labeled with its spectral type and errors 
bars representing the uncertainty of the photometry. The colored symbols are data for the same objects as shown in (a). 
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Fig. 6. — ii's-[3.6] versus [3.6]-[4.5] color-color diagram for all of the sources in our M, L, and T dwarf sample with iC-band photometry, 
excluding known binaries and spectrally peculiar objects. All near-IR photometry has been converted to the 2MASS JHK system using 
the relations of Stephens &; Leggett (2005) (see additional discussion in the text). The plot symbols are the same as those used in Figure 3. 
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Fig. 7. — Similar to Figure 6 except in this case we show i<'s-[4.5] versus [3.6]-[4.5] color-color diagram for the same objects. 
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Fig. 8. — versus _R'-[4.5] color-magnitude diagram for all of the sources in our M, L, and T dwarf sample with i^-band photometry 

and trigonometric parallax, excluding known binaries and spectrally peculiar objects. As in Figure 6, all near-IR photometry has been 
converted to the 2MASS JHK system. Parallax references can be found in Table 1. The plot symbols are the same as those used in Figure 
3. 
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Fig. 9. — Similar to Figure 8 except in this case we show M^g versus J-[4.5] color-magnitude diagram for the same objects. The three 
T dwarfs in the region between the locus of the L dwarfs and the parallel but fainter locus of T dwarfs below have spectral types of TO. 5, 
Tl.O, and T2.0. The sequence of T dwarfs below run from left to right from T3.5 to T8. 
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Fig. 10. — A summary of the trending seen in various color indices versus spectral type in the same spirit as near- and mid-IR plots 
shown in Leggett et al. (2002), Knapp et al. (2004), and Golimowski et al. (2004a). The spectral types are the same as those used in 
Table 1, with optical types being used for the M and L dwarfs and infrared types for the T dwarfs. All J- and JC-band photometry are 
on the 2MASS system, with all transformations made using the relations of Stephen & Leggett (2005). Representative error bars for the 
median errors in the photometry that make up each color index and an error of +/- one spectral subclass are shown in each plot. 



24 



Patten et al. 



11 



12 



^ 13 



14 



15 



16 



(1) 
(a) 

(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 

(10) 2MA0727+1710 

(11) GJ570D (T7.5) 

(12) 2MA0415-0935 (T8.0) 



SDS0151+1244 
SDS0837-0000 
SDS1254-0122 
SDS1750+1759 
SDS0207+0000 
2MA2356-1553 
2MA1546-3325 
2MA0243-2453 
SDS1346-0031 



(TO. 5 
(Tl.O) 
(T2.0) 
(T3.5) 
(T4.5) 
(T5.0) 
(T5.5) 
(T6.0) 
(T6.5) 
(T7.0) 




^ A ' 




m 



m 







[3.6]-[4.5] 



Fig. 11. — A color-magnitude diagram of M3.6 versus the [3.6]-[4.5] color for the T dwarfs in our program with trigonometric parallaxes. 
As in the previous figures in this paper, known binaries and spectrally peculiar objects have been excluded. Each T dwarf plotted is 
numbered 1 to 12 in order of increasing spectral subtype. The accompanying error bars include the estimated photometric errors and the 
published errors in the distance estimates for each individual object (in some cases the error bars are smaller than the plot symbol). The 
blue squares are models for T^fj from 1300 to 700 K (left to right, in steps of 100 K) at g= 10^'^ cm s~^, the red triangles and green circles 
are are models for the same temperature range from at g= 10^ cm s^^ and g= lO*'^ cm s~^, respectively. The corresponding masses range 
from ~15 Mj„piter to 70 M-jupite.r (Burrows et al. 1997). 
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Fig. 12. — A color-magnitude diagram of Ms.g versus the [5.8]-[8.0] color. Otherwise, the format is the same as in Figure 11. The 
theoretical models are also the same and span Teff from 1300 to 700 K, in steps of 100 K, for gravities of 10^'^ cm s^-^ (blue squares), 10^ 
cm s~^ (red triangles), and lO**'^ cm s"'^ (green circles). 
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Fig. 13. — A color-magnitude diagram of M4,5 versus the [4.5]-[5.8] color. Otherwise, the format is the same as in Figure 11. For this 
color-magnitude diagram, the theoretical fits (the same as in Figure 11) are unacceptable. Because CO has a strong spectral feature at 
~4.67 /jm , one can interpret the discrepancy between theory and the IRAC data as an indication that equiUbrium CH4/CO chemistry 
underestimates the abundance of CO in T dwarf atmospheres (GoUmowski et al. 2004a). The magnitude of the discrepancy translates into 
an overestimate in the ~4— 5 /imflux by factors of 1.5 to 3.0. 
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Fig. 14. — A comparison between a theoretical T dwarf spectral model with Tcff /g = 750 K/10^ cm and the IRAC fluxes measured 
for the T7.5 dwarf GJ 570D. A distance of 5.89 parsecs has been assumed (Dahn et al. 2002). The black triangles depict the fluxes 
(with their uncertainties) derived from the IRAC photometry; the horizontal lines indicate the widths of the IRAC bands. This model 
was generated for GJ 570D in 2002 to flt its optical spectrum shortward of 1.0 /^m (Burrows et al. 2002). No further attempt was made 
to improve the flt. A slight (~40%) discrepancy in the 4. 5-/xm bandpass, attributable to a CO abundance excess in the atmosphere, is 
manifest. Despite this, this plot and Figure 7 in Burrows et al. (2002), represent a good fit (that can nevertheless be improved) from 0.6 
/imto ~8.0 /imand indicate that the IRAC data were successfully anticipated. 
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* Unless otherwise noted, optical spectral types are used for M and L dwarfs and infrared spectral types are used for the T dwarfs. 
^ Known binary. 
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TABLE 2 
IRAC Photometry Calibration 



Item 




IRAC Channel 
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3 (b.Sfim) 
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BCD calibration [(MJy/sr)/(DN/s)] 


0.1088 


0.1388 
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0.2021 


Point source calibration [juJy/(DN/s)] 
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7.070 


Aperture correction(error) for 4-pixel radius [mag] 


-0.084 (0.009) 


-0.089 (0.006) 


-0.072 (0.015) 


-0.077 (0.034) 


F.(Vega) [Jy] 


280.9 


179.7 


115.0 


64.13 
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IRAC Photometry and Colors of latb-M, L, and T Dwarfs 
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12 





39 


GJ570D 


T7.5 


13 


80 


(0 


04) 


5 


12 


12 


(0 


02) 


5 


12 


77 


(0 


11) 


5 


11 


97 


(0 


07) 


5 


1 


68 


-0 


64 





80 


2MA0415-0935 


T8.0 


14 


10 


(0 


03) 


5 


12 


29 


(0 


02) 


5 


12 


87 


(0 


07) 


5 


12 


11 


(0 


05) 


5 


1 


82 


-0 


58 





76 


e Ind BC 


T1.0+T6.0 


9 


97 


(0 


01) 


10 


9 


44 


(0 


02) 


10 


9 


39 


(0 


03) 


10 


8 


98 


(0 


04) 


10 





53 





04 





41 



Note. — (1) Source saturated in all four IRAC channels in 30-second FRAMETIME data, 2-second FRAMETIME data used instead, (2) GTO program FID 33 target, 2-second FRAMETIME 
data, (3) GTO program FID 33 target, 0.6-second FRAMETIME data, (4) Channel 1 saturated for 30-second FRAMETIME data, used 2-socond FRAMETIME data for this channel, (5) Channel 
2 saturated for 30-second FRAMETIME data, used 2-second FRAMETIME data for this channel, (6) GTO program PID 33 target, (7) Target source aperture possibly contaminated by flux from 
another nearby source, (8) Target in crowded field, some contamination of source aperture by other nearby sources possible, (9) Spitzer AOR target name incorrectly reads "2MA1410-t-1804", 
(10) Wings of GJlOOl PSF may contaminate source aperture for target, however 2-second FRAMETIME data for channels 1 and 2 agree well with the 30-second FRAMETIME data, (11) Strong 
nebulosity in background, especially in channel 4. 

^ Known binary. 



